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We propose and theoretically investigate spin superconducting qubits. Spin superconducting qubit 
consists of a single spin confined in a Josephson junction. We show that owing to spin-orbit interac- 
tion, superconducting difi'erence across the junction can polarize this spin. We demonstrate that this 
enables single qubit operations and more complicated quantum gates, where spins of different qubits 
interact via a mutual inductance of superconducting loop where the junctions are embedded. Recent 
experimental realizations of Josephson junctions made of semiconductor quantum dots in contact 
with superconducting leads have shown that the number of electrons in the quantum dot can be 
tuned by a gate voltage. Spin superconducting qubit is realized when the number of electrons is odd. 
We discuss the qubit properties at phenomenological level. We present a microscopic theory that 
enables us to make accurate estimations of the qubit parameters by evaluating the spin-dependent 
Josephson energy in the framework of fourth-order perturbation theory. 

PACS numbers: 



I. INTRODUCTION 

Potential benefits of a quantum computer: secure com- 
munication, fast database searching, efficient prime fac- 
torization P, [2, 3 have inspired significant research ef- 
forts. Building a quantum computer requires the real- 
ization of qubits as its elementary units. Useful qubits 
should satisfy two conditions: they can be manipulated 
and read before the quantum information stored in their 
state is lost, and they allow for engineering of a control- 
lable interaction between them. Designing and realizing 
such qubits defines the focus of most research in the field. 

Since the electron spin provides the simplest exam- 
ple of a coherent two-level system, that is, a qubit, 
spin-based qubits very soon became a subject of intense 
theoretical [j| and experimental research. Experiments 
proved relatively long Ti and T2 times for single elec- 
tron s pin s trapped in quantum dots 0, H, 0, , dia- 
mond [l^, [HI, [13, [3 and other materials. Good isolation 
from the environment protects from decohcrence at cost 
of hampering qubit control and read-out. While cur- 
rent research successfully addresses these shortcomings 
[1^ . [isl . [igI . [TtI , the realization of controllable interaction 
between pairs of spin qubits has been so far obstructed 
by numerous practical problems [31 • 

Superconducting qubits do better in this respect. 
Quantum logic gates involving the controllable interac- 
tion of two qubits have been demonstrated in a variety 
of setups [1^ [13, [m, [H, [2^ . Superconducting qubits ex- 
ploit Josephson effect and Coulomb blockade, both im- 
mediately related to electric variables of flux and charge. 
This allows for easy integration of qubits into electric 
circuits and is the reason of the better performance, sen- 
sitivity to external noise that results in relatively short 
decoherence times [22l. [2^. [25| . 

In this article we propose a design of superconducting 
spin qubit and discuss its feasibility and advantages. A 
spin is trapped in a quantum dot connected to super- 
conducting leads, thereby forming a Josephson junction. 



see Figure [TJ Owing to spin-orbit interaction, the super- 
conducting phase difference polarizes the spin. This pro- 
vides means to read and manipulate its quantum state. 
We detail the operation of a single qubit and the design 
of qubit-qubit interaction that allows to make quantum 
logic gates. In particular, we emphasize the prospect of 
all electrical manipulation of the qubit state and qubit- 
qubit coupling. To prove feasibility of the design, we 
present microscopic calculations and numerical estimates 
of the spin- and gate voltage-dependent Josephson en- 
ergy. 




FIG. 1: Spin superconducting qubit. a) Quantum dot with 
odd number of electrons is connected to superconducting leads 
(L and R) biased at phase difi'erence b) Energy levels in 
quantum dot and the leads. 

Quantum dot connected to superconducting leads is an 
essential element of our design. Superconducting quan- 
tum dots have received theoretical attention rather early 
[26l . [29j . About ten years later, experimental break- 
through has been achieved by making good contacts 
between semiconducting nanowires and superconducting 
leads [13, [m . Gate electrodes put close to the nanowire 
can be used to create potential barriers in the nanowire, 
thereby defining a quantum dot. The quantum dot is 
in the Coulomb blockade regime, that is, the number 
of electrons is tunable. It has been proven that such 
a dot can be included into a superconducting circuit as 
a Josephson junction carrying a supercurrent [3l|. The 
idea of our proposal is to use this setup by keeping odd 
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number of electrons in the dot. In this case, the ground 
state of the resulting Josephson junction is a spin dou- 
blet. This is advantageous in comparison to an earlier 
proposal concerning Andreev quantum dot [s^ . where 
the spin doublet corresponds to an excited state of the 
system. Spin-orbit effects in quantum constrictions and 
dots have been discussed in [23, [2^ . 

The structure of the article is as follows. In section 
ITTl we describe the qubit phenomenologically, discussing 
single qubit manipulation and design of qubit-qubit inter- 
action. Section|TTT]focuscs on the microscopic description 
of the Josephson junction with a single spin, leading to 
expressions of the spin-dependent and spin-independent 
parts of the Josephson energy in terms of the junction 
parameters. Section IIVI discusses different parameter 
regimes, providing order of magnitude estimations of the 
Josephson energy. Section |V] presents numerical calcula- 
tions of the relevant quantities. Section IVTl concludes. 



II. PHENOMENOLOGICAL DESCRIPTION 

We proceed with the phenomenological description of 
the qubit, postponing for now microscopic analysis. An 
important feature is that the superconducting current 
flowing through the Josephson junction is determined by 
the state of the spin enclosed. This can be understood 
as follows. Junction current is the result of transfer of 
Cooper pairs between the leads via coherent tunneling 
events. Conventionally, tunneling events are spin con- 
serving. However, strong spin-orbit coupling mixes the 
spin states of electrons as they tunnel between orbitals 
in the superconductor and orbitals in the dot [s^ . The 
resulting non-spin-conserving tunneling amplitudes de- 
pend strongly on the wavefunction of the levels involved. 
When a Cooper pair tunnels via two different levels in the 
dot, the initial spin-singlct configuration acquires com- 
plex spin structure, due to the different spin-dependent 
tunneling amplitudes. If one of the two levels involved is 
occupied by a single electron, its spin restricts the path- 
ways of Cooper pair tunneling, due to Pauli exclusion. 
This directly couples current to spin, lifting the degener- 
acy between spin up and spin down states. 



A. Effective Hamiltonian 

We describe the spin polarization effect by the follow- 
ing junction Hamiltonian: 



H = 



Ej cos [if) + {Eso 



sin {lp) 



(1) 



where a is the spin operator. The pseudovector e*so de- 
fines the polarization axis in three dimensions. Its mag- 
nitude and direction depends on the spin-dependent tun- 
neling amplitudes as well as on the positions of the levels 
in the quantum dot. As a result, Eso is independent of the 
superconducting phase difference Lp, but does depend on 



gate voltage, as the gate electric field modifies quantum 
dot wavefunctions and levels. 

We will present estimations of Ej and e^o in scction lTVl 
For present purposes it is enough to assume that typical 
values of eso are somewhat smaller than Ej . Actually, Ej 
can be made zero by a certain choice of gate voltages. 
In the vicinity of this point, e^so may be bigger than Ej. 
However, we do not concentrate on this case. 

The junction forms a qubit: it may be found in two spin 
states. It) and ||), that differ in spin projection along the 
polarization axis eso- At fixed phase, the energy spacing: 
/S.E = 2 I ego I I sin ((^)|, leads to a measurable difference in 
superconducting current (Figure [2]): 



A/ = ±{2e/h) |cos(^)| 



(2) 



The sign is determined by the direction of the spin along 




FIG. 2: Two states of a simple superconducting circuit (Fig- 
ure[3K) containing the qubit junction, a) Energy as a function 
of phase drop over the qubit junction Lp. b) Loop current /(i/j) . 



This provides the means to read the qubit state. 




FIG. 3: Two simple superconducting circuits, a) Supercon- 
ducting loop interrupted by qubit junction (characterized by 
Eji and e*Bo) and conventional Josephson junction (character- 
ized by Ej2)- The phase of the two junctions are modulated 
by the magnetic flux: ^1 4- = 2-K<f)\oop/ b) Two super- 
conducting loops are presented, each interrupted by a qubit 
junction. The loops have a common side that is interrupted by 
a conventional Josephson junction Ej^. Magnetic flux flowing 
trough the loops is represented by (^loop 1 and (/)ioop 2- 

To illustrate, consider a simple superconducting circuit 
consisting of a loop interrupted by the qubit junction 
and a conventional Josephson junction, see Fig. [3^. We 
apply magnetic flux through the loop (pxoop- When the 
Josephson energy of the conventional junction is much 
larger than the Josephson energy of the qubit Ej2 3> Sji, 
the phase induced by the magnetic flux is acquired mainly 
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by the qubit junction Lp2 ~ 0. The phase difference over 
the qubit junction can be fixed by fixing the magnetic 
flux Lpi w 27r0ioop/'?!>Oi 00 = 7rft/e being the magnetic 
flux quantum. The current through the loop is plotted 
as a function of phase for each of the two states of the 
qubit (see Figure[2]). The magnetic flux generated by this 
current can be measured by a nearby SQUID loop, not 
shown in Figure [31 This is a common technique used to 
measure the state of superconducting qubits [23| ■ 

The maximum value of the induced flux is achieved 



when E- 



Ej2 ~ |e*so| and is of the order of the mag- 
netic flux quantum (f>Q. It is interesting to compare the 
induced flux, with the flux generated by the magnetic 
dipole moment of the electron spin confined in the quan- 
tum dot. Considering an area element perpendicular to 
the dipole moment and situated at the distance r from 
it, the flux flowing through the area scales with the dis- 
tance as 4> oc . We find the distance rg where the 
magnetic field of a single electron produces a quantum of 



magnetic flux via the relation e /(47reor'o) 



This 



distance is well into the subatomic region, in the order of 
ro « 10-1^ m. 



B. Spontaneous Currents 

Let us consider another regime, which is useless for 
qubit applications, but interesting from the point of view 
of general physics. Let us consider vanishing flux through 
the loop and assume formula Eji,Ej2 <SC |e*so|- In this 
case, the phase dependent energy Eq. ([T]) can be ex- 
panded at small (p 

E{^)^iE,,+E,2)^V2±\e,,\^ , 

with opposite signs corresponding to opposite spin orien- 
tations. 




b) 




FIG. 4: Superconducting loop without magnetic flux. The 
dependence of energy on phase shows that the minimum en- 
ergy is achieved at non- vanishing phase ±(fio- The sign of the 
equilibrium phase depends on the orientation of the spin in 
the qubit junction. 



zero flux applied, the system is time-reversible. The sit- 
uation just described can be envisaged as spontaneous 
breaking of time reversibility. Indeed, the energies of the 
states with opposite spins are precisely the same. We will 
present detailed description of the situation elsewhere. 



C. Single Qubit Manipulation 

Let us turn our attention to manipulation of the qubit 
state. It is common to use pulses of an ac field of res- 
onant frequency w = 2 |e'so| /S |sin((/3)|. Magnetic fields 
perpendicular to the loop plane induce modulations of 
superconducting phase, but do not allow resonant ma- 
nipulation as this does not change ego- 

In-plane magnetic fields can polarize qubit spin, de- 
flecting it from the direction of Cso- However, this ap- 
proach is difficult to realize experimentally due to mis- 
alignment. The magnetic field component perpendicular 
to the loop contributes to the total magnetic flux, chang- 
ing the properties of the qubit at the same time as the 
manipulation is performed. 

Fortunately, magnetic field is not needed. Rabi os- 
cillations can also be induced electrically via the gate 
electrodes. We remind that the direction of Ego depends 
on the position of levels in the dot. Therefore, shift- 
ing the gate voltage would also rotate e*so- To illustrate, 
let us assume that a change of the gate voltage leads 
to a corresponding change in polarization pseudovector 
e*so e*so + i^e'so. We describe the effect of the resonant 
pulse by the time-dependent qubit Hamiltonian: 

H{t) = [\eso\ a^ + Seso-a cos (cjt)] sm{ip) , (3) 

where the z-axis is chosen along e*so- Assuming \Seso\ ^ 
|eso|, the dynamics of the qubit is described by Rabi os- 
cillations with frequency: 



^so,y) I sin {(f) 



We stress that if the change Je'so is parallel to the initial 
direction of e*so Rabi oscillations do not occur. In this 
case the qubit retains its initial state. Thus, change of 
Scso in the perpendicular direction is essential for qubit 
manipulation. 

We also mention that since electric fields are easier to 
localize in space than magnetic fields, electrical manipu- 
lation is also advantageous in view of controlling qubits 
individually. 



D. Design of Qubit-Qubit Interaction 



We see (Figure |4]) that the equihbrium superconduct- 
ing phase is non-zero. It takes opposite values depending 
on spin orientation. The current flowing through the 
junction is also non-zero / = Ej2 (po- We stress that at 



Further, let us focus on qubit-qubit interaction. Two 
qubits can be included in a superconducting circuit, such 
that the spin-dependent supcrcurrents flowing through 
the two qubits interact magnetically. The interaction 
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does not modify the polarization pseudovectors. As a 
result, the qubit-qubit interaction is of the Ising type: 



H = Hidl + H2(J2 + alal 



(4) 



choosing the z-axis along ego for each qubit. The Ising- 
type interaction is sufficient to perform the CNOT opera- 
tion, which, in combination with single qubit operations, 
enables universal quantum computations [36j . 

The CNOT gate is an operation on two qubits which 
has the effect of changing the state of one (target) qubit, 
only when the other (control) qubit is in the excited 
state. We propose a realization of CNOT gate using non- 
oscillating pulses of Hi2/h of length r. The pulse shifts 
the relative phase between two qubit states, as follows: 
states Itt) and |J,|) gain phase factor exp{iHi2T/h), while 
states Iti) and ||t) gain phase factor exp{—iHi2T/h)). 
Tuning the length of the pulse such that Hi2T/h = 7r/4, 
one obtains a phase-shift gate Gphase that can be com- 
bined with Bloch sphere rotations by tt/2 of the control 
and target qubit around the coordinate axes Ra;_y,2(7r/2), 
to achieve the controUed-NOT gate: 



CNOT = e"/4 Rii)(7r/2)R[2)(-7r/2) 

Ri2)(7r/2)GphascR|,'n'r/2) , 



(5) 



where the superscripts (1) and (2) imply that the oper- 
ation is performed on the control and target qubit, re- 
spectively. 

We have the opportunity to tune the interaction and 
also to switch it on and off. To give an example, consider 
the following circuit: two qubits with their superconduct- 
ing loops connected in parallel (Fig. [5)3). The interaction 
energy can be tuned by the following parameters: the 
magnetic flux flowing through the superconducting loops 
('/'loop 1 and 0ioop 2) and the Josephson energy of the con- 
ventional Josephson junction Ej^. To lowest order in the 
polarization pseudovectors |frsoi.2|, the interaction takes 
the form: 



H12 



|esoi| |eso2| 



cos {(p^ + ipi) COS {ip^ + ip2) ,(6) 



where the complex-valued energy E denotes: E = 
Ejie^'^^+Ej2e^'^^+Ej3, and ip^ denotes its complex argu- 



ment 'fi^= tan 



The angles (pi,2 are pro- 



ImE/ReE 

portional to the magnetic fluxes: Lpi^2 = 27r (0ioop i,2/0o), 
with (j)Q the magnetic flux quantum. Setting <pi = </?2 = 
7r/2, it is possible to tune the interaction simply by 
controlling Ej^. For Ej^ 3> {Eji + Ej2) the interac- 
tion is turned off H12 0, while in the opposite limit 
Ej3 <^ {Eji -\- Ej2) the interaction achieves a maximum 
H12 = — Icsoil |eso2|/|-^|. Such control of the conven- 
tional Josephson energy can be achieved u sing exclusively 
electrical means, as was demonstrated in 1311. 



III. MICROSCOPIC DESCRIPTION 

A. Hamiltonian 

The total Hamiltonian comprises terms describing 
leads (left and right), quantum dot and tunneling be- 
tween electronic states in the leads and in the dot: 



QD 



(7) 



The electronic states of the left lead are labeled by / and 
spin index a and are affected by the superconducting 
order parameter Ae*'''^ . 



aia 



Ae^^-g^-'alal+g'^'^'Ae-^^-ai^a,,, , (8) 

were (a) denotes the creation (annihilation) operators 
and ^ denotes the energy of the levels in the normal metal 
state, counted from the chemical potential of the lead. 
The matrix g = ia^ ensures that the two electrons are 
a spin singlet when forming a Cooper pair in the leads. 
In the superconducting state, the energy of the levels is 

Similarly, electronic states in the right lead are labeled 
by r and are affected by Ae*"^^ . 

The Hamiltonian of the quantum dot contains charg- 
ing energy term along with the terms describing non- 
interacting electrons in levels labeled by m, ^m0.l„a,n, 



Ec 



(9) 



with N being the operator of the total number of elec- 



trons N 



al^ttm. No = VgCg/e represents the effect 



of gate voltage which allows to tune the total number 
of electrons in the dot. Cg here is the capacitance to 
the gate. Since it only appears in combination with gate 
voltage, we find it convenient to rescale the gate voltage 
Vg Cg/C Vg, wlth C bclug the total capacitance of the 
junction. In these units, Nq = eVg/Ec- 

We stress that all these terms are not affected by spin- 
orbit coupling. The reason for this is that spin-orbit 
coupling retains the double degeneracy of the electronic 
states. Instead of spin doublets they become Kramers 
doublets, so index a in this case refers to the compo- 
nents of Kramers doublet, rather than to the original spin 
state. The only place where spin-orbit coupling plays a 
role is the tunneling part of the Hamiltonian and can be 
included into spin-dependent tunneling amplitudes (ssj : 

= ^ TZal^aip + TZaLa^np + H.c. ,(10) 

/,r,m,p,fT 
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where p and a are spin indexes. 

As a result of time reversal symmetry, we may write the 
tunneling matrices in spin space: 



TZ = rO„<5,.+z ^tL<- (11) 



are chosen to be real. The spin- 



where and t]^ 
independent part is symmetric Tj'^^ = and spin- 
dependent part is antisymmetric = — with respect 
to interchanging the state indices. 

In principle the tunneling matrix elements can be eval- 
uated by computing overlap of spin-dependent wavefunc- 
tions of corresponding components of Kramers doublets. 
It is known that wavefunctions become random either 
because of scattering at defects in the leads and in the 
dot, or because of chaoticity of classical electron motion 
in the dot. If we consider extended system, the statistics 
of spin-dependent part of such overlaps is governed by 
length scale Zgf that is a spin flip length, induced by spin- 
orbit scattering. At distances exceeding this scale, spin 
orientation of Kramers doublets becomes completely ran- 
dom, so that r° is of the order of |t| . The length scale Zgf 
corresponds to time scale Tgf , spin-flip time. If we apply 
this now to a quantum dot, the spin structure becomes 
random for each level in the dot, provided the mean level 
spacing in the dot 5s is comparable or smaller than h/Tsi. 
The ratio of the two defines the strength of spin-orbit 
coupling. For our estimations, we assume t/T'^ « 0.1 
and random distribution of the direction. Preliminary 
experiments with InAs nanowires [37j confirm this by or- 
der of magnitude. 



B. Josephson Energy 

In the following we provide a microscopic description 
of the Josephson effect in the presence of spin-orbit cou- 
pling. We treat Cooper pair transport through a quan- 
tum dot in the regime of Coulomb blockade, with no 
bias voltage and disregarding effects of thermal excita- 
tion. Within the Coulomb diamonds, the only mecha- 
nism of transport is cotunneling. Four tunneling events 
are sufficient to transfer a Cooper pair between the leads 
[3^ and these processes dominate for sufficiently small 
normal state conductances Gn ^ Gq = /nh. We em- 
ploy fourth-order perturbation theory in the tunneling 
amplitudes and find the Josephson energy as the pertur- 
bation correction to the ground state energy. 

The ground state of the unperturbed Hamiltonian 
Ha = i?L + ^^R + ^^QD is degenerate with respect to the 
number of Cooper pairs on the superconducting leads. 
The wave function can be written as a direct product of 
three wave functions: 1/' (A^l, A^r) = V'l (A^l) ® V'QD ® 
V'R (-^l) corresponding to the two leads and the quan- 
tum dot. A'l and A^r represent the number of electrons 
in the respective leads. 



In the presence of tunneling, the degeneracy of the 
ground state is lifted and the energy splitting between 
states with different number of Cooper pairs in the leads 
corresponds to the Josephson energy: 



V'(-2,2) 
V.(2,-2) 
2Re{(v^( 



f(4) 

-2,2^ 



V'(0,0)^e*'^ + 
V^(0,0)\e-^'^ 



(4) 



V'(0,0))e^'^}(12) 



where A^l and A^r have been set to zero. 

Here, operator T*^^) represents the correction to the 
amplitude of the fourth order in tunneling Hamiltonian: 



f(4) =i7T- 



1 



1 



1 



£^0 — Hn Eq — Ho Eq — Ho 



Ht ■ (13) 



The structure of the fourth-order correction is as fol- 
lows. Each tunneling operator Ht describes the tran- 
sition of one electron between a state in the dot and a 
state in one of the superconducting leads. There are 24 
distinct sequences resulting from the permutations of the 
four tunneling events considered. Each process involves 
three intermediary states. The energies of the virtual 
states appear as the three denominators of Eq. [T31 The 
three energy denominators unambiguously characterizes 
the sequence of tunneling events. 

In contrast, the spin structure does not depend on the 
order of individual tunneling events. The spin-structure 
of T^^^ can be recovered from the product of four tunnel- 
ing amplitudes describing hopping of the two electrons 
between the leads and the quantum dot: 

l,r 

'^^TmiTinTnrTrm 5 (^L — 6) ^ i^B, — £,r) = 



l.r 



, (14) 



where m and n label states of the quantum dot. 

Let us explain the properties of the two terms in Eq. 
(|14p . In the following, we neglect terms proportional 
to the square of the small term t/T^ . Since the leads 
have continuous spectra, it is convenient to introduce 
transport rates which are continuous functions of en- 

2 



5{^-£.Kr))- The 



ergy: TL{R).m{i) = ^Ti/h 2J ['^m/W 

l(r) 

spin-independent term {^l,£,r), has different prop- 
erties in the case m = compared to the case m 7^ n. 
Terms diagonal in the dot levels have the following sim- 
pie form: P°^(a,?fl) = {h/2TrfT L,m{iL)T R^mi^R). As 
a result, they are always positive. In contrast, the sign of 
Prnn {S,l,£,r) when m n can be related to the parities 
of the wavefunctions of the two dot states. In analogy to 
the case of a particle in a box, the parity of wavefunctions 
in a quantum dot alternates between neighboring states. 
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As a result, the sign of n^m (CijCfi) alternates and 
the sum over all states in the quantum dot averages out: 

The contribution corresponding to the spin-dependent 
term in Eq. (|14|) changes sign for different spin orien- 
tations. Therefore, this term cancels out if the levels m 
and n are both cither empty or filled with electrons, as 
there is no asymmetry between spin up and spin down 



terms. However, if either state m or n is the state filled 
with a single electron, the contribution is finite, giving 
rise to the spin polarization effect. 

After calculating the energy denominators for each 
tunneling sequence and integrating over all states in the 
leads and all pairs of states in the dot, we obtain the fol- 
lowing result for the spin-independent Josephson energy 




E, = 



A2 



Ml, A/2 



2 E ^Mrn {^L,^r) Jeh{^L,£,R, ^M, £,7,1) + ^ ^MO (^ij Ci?,) {Jee {£l,£,R, £,M , £,o) - Jeh {£l,£,R, £,M,£o)) 



M.m 



M 



E ^Om (C-Li ^fl) (Jhh {S,L,£r, ^Oi Cm) " Jeh {^Li^R, Cm)) " 2Pgo (Ci, C-R) Jeh{^Li£,Ri C07 Co) 



(15) 



and for the polarization pseudovector: 



A2 r d^ 



PmO {^L,£,r) {Jee {^L , £,R, £,M , £,q) + Jeh {^L , £,R, £,M , £,o)) " 



M 



E {£.L,£.r) {Jhh {£,L,£,R,^0,£,m) + Jeh {£,L, £,R, £.0, £.m)) 



(16) 



where to, mi and TO2 label filled dot states, M, Mi, 
M2 empty dot states and labels the half-filled state. 

We have used the notation e^n = y^A^ + Ci The J- 

functions contain the energy denominators and arc dif- 



ferent for processes where the Cooper pair is transfered 
via two electrons, two holes or one electron and one hole, 
respectively: 
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1 



1 



1 



1 



1 



1 



1 



EL + Cm + E{e) EL +EnEii + + E{e) EL + in + E{e) EL + ER ER + Cm + E{e) 

11 111 



EL + Cm + £^(e) Cm + Cn + E{2e) ER + U + E{e) EL + Cn + E{e) Cm + Cn + E{2e) ER + C,„ + E{e) 

111 111 

eL + Cm + -E(e) Cm + Cn + E{2e) ER + Cn + E{e) ^ EL + Cn + E{e) C,„ + Cn + E{2e) er + Cn + E{e) 

Jhh{iL,£,B., Cm, Cn 



1 



CL - Cm + E{-e) EL +ERER- Cn + E{-e) EL - Cn + E{-e) EL + ER Er - Cm + E{-e) 

11 111 



CL - Cm + S(-e) -Cm - Cn + £^(-26) Er - Cm + E{~e) El - Ca + E{-e) ^Cm - Cn + E{-2e) ER - Cm + E{-e) 

111 111 



EL - Cm + S(-e) -Cm - Cn + £;(-2e) Cfl - Cn + -B(-e) et - Cn + E{~e) -Cm - Cn + E{-2e) ejj - Cn + £^(-e) ' 
Je/.(CL,Ci?,Cm,Cn) = + + E{e) El + Er El - Cn + E{-e) + ei^ + Cm + £;(e) eL + efl - Cn + i?(-e) + 

111 111 



+ Cm + E[e) Cm - Cn + + e_R ER + Cm + E{e) ER - Cn + E{-~e) Cm - Cn + EL + Cfl e_R + Cm + i?(e) 



CL + Cm + £^(e) Cm - Cn + CL + e_R EL - Cn + -B(-e) Er - Cn + -B(-e) Cm - Cn + + Cfl CL - Cn + -E'(-e) 



IV. ESTIMATIONS 

Let us estimate the typical magnitude of Cso and com- 
pare it to the magnitude of spin-independent Josephson 
energy. We will see that the relative magnitudes as well 
as its absolute value cannot be just estimated by the typ- 
ical strength of spin-orbit coupling |tj /Tq- The estima- 
tions depend on three energy scales in the problem: 5s , 
Eq and A. Besides, it depends on the energy distance 
to the diamond edge. 

It has been already shown in [33|] that Josephson en- 
ergy and critical current exhibit spectacular peculiarities 
at the diamond edge. These peculiarities are no singu- 
larities. This is related to the diamond structure in the 
presence of superconductivity [s^l- There is a bistabil- 
ity region in interval of width A(eVg) = 2A around each 
diamond edge. The singularities in denominators of per- 
turbation theory occur at the edges of bistability region 
for the charged state of higher energy. We always as- 
sume the dot to be in the lowest energy state. This 
saves us from singularities. However, the proximity to 
the point of singularity gives a spectacular increase or 
decrease in Josephson energy at the diamond edge. To 
account for this in our estimations wc introduce the min- 
imal energy distance to the edge of Coulomb diamond 
E± = I mm{E{e), E{—e)) \ > A. We envisage two sepa- 
rate situations: E± « Ec and E± ^ Ec- 

The dominant contribution to the Josephson energy 
consists of the terms with smallest energy denominators. 
The magnitude of the energy denominators depends on 
the energy of the dot levels involved in tunneling. We 
find that the energy interval for quantum dot levels in- 
volved in the dominating contribution has the width E± , 



even in the regime E± ^ A. Let us define the number 
of levels within this interval Ns as the integer part of 
E±/Ss + 1. If (5s > E± the Josephson energy is domi- 
nated by the contribution of a single level. If Ss < E± 
there are multiple levels participating in the dominating 
contribution. 

Using Eq. (fT5|) we find the following estimate for the 
spin-independent Josephson energy: 

J, ^ M ^^^^ 

Similarly, we can estimate the spin-dependent component 
using Eq. (fTB)) : 



' ^ ^ A E±{miix{E±,5s)) Tq ' 

In all parameter regimes considered, both Ej as well 
as |eso| increase as E±/A decreases, i.e. as we approach 
the edges of the Coulomb diamond. We can explain this 
in terms of the energy of the intermediary virtual states. 
In the middle of the Coulomb diamond the energy cost of 
adding an electron to the quantum dot is maximum and 
the high energy of the intermediary states reduces the 
probability of Cooper pair tunneling. Toward the edges 
of the diamond, tunneling processes involving the state 
closest to resonance will have intermediary states that 
are lower in energy, resulting in the increase in Josephson 
coupling. 

Let us turn our attention to the multi-level regime, 
where Josephson tunneling is a result of the interference 
of tunneling processes that involve all pairs of the Ns rel- 
evant states. The spin-dependent and spin-independent 
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terms do not scale in the same way with the number 
of levels involved. In the case of Ej, the dominant 
contribution results from summation of diagonal spin- 
independent elements ^'mm(^i"^fl)- ^ result, Ej 
scales with the number of levels Ns- In the case of e'so, 
the result of summation over tunneling contributions is 
equivalent to the distance traveled in a random walk in 
three dimensions, after Ns steps. The average result in 
this case scales as \^Ns- 

The ratio of spin-dependent to spin-independent con- 
tributions is: 



1 



E± 



E, 



/]V^ mSix{E±,6s 



To 



Apart from the trivial dependence on spin-orbit coupling 
strength, the ratio |eso| /Ej is further reduced by factors 
that depend on the gate voltage. We distinguish two im- 
portant regimes: the single level regime 63 > Eq-, and the 
multi- level regime 5s ^ £'c, where multiple dot states 
contribute. 

In the single level regime we find |e'so| /Ej w E±/Ss- 
Thus, it is large in the middle of the diamond, where it 
scales as Ec/Ss, decreasing towards the edge, where it 
scales as A/ 63- The behavior has a simple explanation in 
terms of energy of intermediary virtual states. The dom- 
inant contribution to Ej arises from processes involving 
tunneling of both Cooper pair electrons through the dot 
level closest to resonance. Such processes do not con- 
tribute to the polarization effect; as we have seen, only 
processes involving two different levels contribute. Thus, 
the energy of intermediary states decreases faster as we 
approach the diamond edge for dominant processes con- 
tributing to Ej . 

We can also estimate the dependence of the direction of 
e'so on the gate voltage. For this, we define the change of 
e*so in the perpendicular direction, in the following way: 

This choice is motivated by the desire 



to achieve electrical manipulation of the qubit. We have 
shown in sectionllTlthat change of e'so perpendicular to its 
initial direction is an essential ingredient for observation 
of Rabi oscillations. 

We find the following estimate of the perpendicular 
derivative: 



de*s 



Tr.r 



Li R 



£:± [max(£;±,(5s)]^ To 



To understand why the direction of e'so depends on gate 
voltage, we study the structure of Eq. ([TB)) . We have 
mentioned that e'so is the sum of Ns vectors characteriz- 
ing spin-dependent tunneling via two different states in 
the dot, one being the half-filled state. The weight of 
each vector is determined by the energy denominators, 
and thus depends on gate voltage. It is thus necessary to 
include the contributions of at least three different lev- 
els, one being the half-filled state, in order to estimate 
the change of the orientation of e'so- As a result, for large 



level spacing the derivative in the perpendicular direction 
is reduced by a factor of {E±/5s)^ ■ In comparison, the 
derivative of the modulus of e'so decreases slower for large 
level spacing, only as E±/Ss- 



V. NUMERICAL ANALYSIS 

We perform numerical calculations of the spin- 
independent and spin-dependent parts of the Josephson 
energy, based on Eq. (fTS)) and p6|) . We are particularly 
interested in the ratio | e'so | / Ej , relevant for experiments 
aiming to measure the separation between qubit states, 
and on the dependence of the direction of e'so on gate volt- 
age, relevant for experiments aiming to perform electrical 
manipulation of the qubit. 

To set up the calculations, we assume that the tunnel- 
ing amplitudes are independent of the lead states. We ac- 
count for the fact that experimental realization of quan- 
tum dots docs not permit control of the resulting local- 
ized states: we choose random energy spacing between 
the dot states and random absolute value of tunneling 
amplitudes. The parities of the localized wavcfunctions 
are also chosen randomly, resulting in a random sign as- 
sociated to off-diagonal spin-independent contributions 
Tm n^m {^l,£r)- The parameters characterizing the dot 
levels are the average level spacing Ss and average modu- 
lus of the spin-independent and spin-dependent tunneling 
amplitudes, respectively (|T°|) and As estimated, 

the ratio |eso| /Ej is proportional to and we 

include this ratio in the energy unit of the spin-dependent 
term. 

Additional parameters in the calculation are the su- 
perconducting energy gap 2A and the charging energy 
Ec- For the numerical analysis, we need to consider a 
finite number of levels of the quantum dot. The results 
presented are obtained including a number of = 20 
quantum dot levels. 

We vary the gate voltage over a large domain, per- 
mitting observation of multiple Coulomb diamonds (see 
Figures [5][5]) . The size of the diamonds observed is Ec in 
the case of odd number of electrons and increases by the 
level spacing in the case of diamonds with even number 
of electrons. 

The units in Figures [5][5] are chosen in accordance with 
the estimations presented in the section above, such that 
the value of unity corresponds to the estimated value of 
the quantity close to the edge, i.e. at E± — A. The 
results confirm our estimations. 

We focus on four regimes differentiated by high charg- 
ing effects Ec/ A. = 10, see Figures [5] and [6l and relatively 
low charging effects Ec/ A = 1.5, see Figures [7] and [51 We 
also compare large average level spacing 5s /Ec = 1-5, 
see Figures [S] and [71 to regimes where the level spacing is 
smaller 5s /Ec =0.1, see Figures [6] and [8] 

A general feature of the results in Figures [M51 is that 
both spin-independent and spin-dependent terms in the 
Josephson energy increase as we approach edges of the 
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FIG. 5: Numerical results obtained for the regime: Ec I ^ = 
10, 5s /Ec = 1.5. Six diamonds are presented: diamonds 
with even number of electrons are represented by lighter shad- 
ing, while darker shaded regions represent diamonds with odd 
number of electrons. For details regarding the quantities plot- 
ted, see Section FVl 
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FIG. 7: The same as in FigureO for the regime: Ec /A = 1.5, 
5s/Ec = 1.5 
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FIG. 8: The same as in Figure[5l for the regime: Ec /A = 1.5, 
5s /Ec = 0.1 



FIG. 6: The same as in Figure[5l for the regime: Ec/ ^ = 10, 
5s /Ec = 0.1. 



diamond. Furthermore, comparing Figure [7^, where 
Ec ^ A, with Figure [5]i, where Eq ^ A, we can con- 
clude that for lower charging energies the ratio between 
Josephson energy in the middle of diamond and energy at 
the edge is reduced, in agreement with our estimations. 
In Figures [THHK Ej is shown to decrease in modulus to- 
wards the edge. This an unusual behavior that can be 
explained if there is a point in the higher energy state of 
the bistable region E± < A where Ej would change sign. 

Let us focus on the spin-independent contribution to 
the Josephson energy, presented in panel A of Figures [5]- 
[51 In the single level regime, we expect Ej to be positive 
for diamonds with even number of electrons and nega- 
tive for odd occupancy, as explained in [3l|. This be- 
havior, observed in Figure [5^ in the first two diamonds, 
would dominate for Ss/Ec ^1. In the regimes pre- 
sented, an interesting feature occurs: there are regions 
of gate voltage where Ej changes sign within a single 



diamond (point of supercurrent reversal [3ll|). As a re- 
sult, Ej is dramatically suppressed. In these regions, 
I ego I may provide the dominating contribution to the total 
Josephson energy. This can be observed in panel D when 

|eso| I \J E'j + Icsol^ ~ 1, see Figures [5] and [6l Comparing 
the different regimes, we observe that regions where Ej 
changes sign are more likely to appear if 5s ~ Ec- In 
addition, the larger variations of the Josephson energy 
observed for Ec ^ A compared to the regime Ec ^ A, 
further increase the probability of sign reversal. 

Turning to the polarization vector ego, panels B and C 
of Figures [5][5] present its modulus and the variation of its 
orientation as a function of gate voltage. The modulus 
|eso| exhibits similar behavior as the absolute value of the 
spin-independent term Ej . In contrast to Ej , we do not 
observe regions of gate voltage where jesol vanishes. It is 
not surprising: the probability that all three components 
of the polarization vector would vanish at the same value 
of gate voltage is very small. 

In panel D of Figures [S][5] we plot the normalized ratio 
between the spin-dependent and spin-independent terms 
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|eso| / Y -E"! + |e*so|^. We note that the ratio mcreases m 
the vicinity of regions where Ej vanishes, as is the case in 
Figures [7][8li. These regions arc exceptional. The com- 
mon case is represented by diamonds such as diamond 1 
in Figure [SJi, where one observes a maximum of the ra- 
tio in the middle of the Coulomb diamond and decrease 
toward the edges. This is in agreement with our estima- 
tions, reflecting that Ej diverges faster than |e*so| as we 
approach the diamond edge. 

Let us focus on the direction of ego and its dependence 
on gate voltage. We study the derivative deso/dVg, pro- 
jected in the perpendicular direction to e'so- As a gen- 
eral feature, the derivative is significantly reduced in the 
middle of the diamonds, as compared to the edge. It is 
interesting to note that the change in the perpendicular 
projection is not as smooth as the variation of eso with 
gate voltage. In the regime A 6s, see Figures [5][7l:, 
we find regions of gate voltage where the perpendicular 
derivative either vanishes, or abruptly changes its behav- 
ior. These features can be explained as follows. In this 
regime Ns ~ 1 near the edges, meaning that the direc- 
tion of is dominated by contributions of few levels. 
It is possible to find regions of gate voltage where the 
weight of the dominating contribution vanishes, similar 
to the case already discussed for Ej. In the vicinity of 
such points, the weight of the dominating contribution 
changes sign. This is represented by the sharp turning 
points observed in Figures [5][71;. In comparison, the case 
Ss presented in Figure[8t shows that the behavior is 
smooth. In this regime the contribution is dominated by 



terms of multiple levels, reducing the probability to en- 
counter values of the gate voltage where the dominating 
contribution to the perpendicular derivative vanishes. 

We may also conclude that ego varies mainly laterally 
when the gate voltage is set to the middle of the diamond. 
Here, the variation of jfTsol vanishes. As the gate voltage 
is tuned towards the diamond edges, the lateral variation 
is overcome by the faster divergence of |eso|- 

VI. CONCLUSION 

In conclusion, we have outlined our proposal of spin 
superconducting qubit. Such unit would combine the 
natural representation of two level system in terms of 
electron spin and advantages of superconducting qubits. 
Spin and superconducting qubits can be operated within 
the circuit, the fiux and spin degrees of freedom can be 
easily entangled. We have demonstrated feasibility of 
all electric manipulation of superconducting qubits and 
more complicated quantum gates made of such qubits. 

The microscopic analysis presented shows that the 
spin-dependent part of the Josephson energy can be made 
sufficiently large, at least for semiconducting devices 
where spin-orbit interaction is intrinsically strong. We 
predict spontaneous breaking of time-reversal symmetry 
in the loops containing spin superconducting qubits. 

We acknowledge fruitful discussions with L. P. 
Kouwenhoven, H. Keijzers, and S. Frolov. This work is 
part of the research program of the Stichting FOM. 
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